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Abstract 

The Chapman-Enskog method of solution of kinetic equations, such as the Bohzmann equation, is based on an expansion 
in gradients of the deviations of the hydrodynamic fields from a uniform reference state (e.g., local equilibrium). This paper 
presents an extension of the method so as to allow for expansions about arbitrary, far-from equilibrium reference states. The 
primary result is a set of hydrodynamic equations for studying variations from the arbitrary reference state which, unlike the 
usual Navier-Stokes hydrodynamics, does not restrict the reference state in any way. The method is illustrated by application 
to a sheared granular gas which cannot be studied using the usual Navier-Stokes hydrodynamics. 



PACS numbers: 05.20.Dd, 45.70.Mg, 51.10.-(-y 



I. INTRODUCTION 



The determination of the one-body distribution function, which gives the probabihty of finding a particle at some 
given position, with a given velocity at a given time, is one of the central problems in nonequilibrium statistical 
mechanics. Its time-evolution is in many cases well described by approximate kinetic equations such as the Boltzmann 
equation^, for low-density gases and the revised Enskog equation0,|3], for denser hard-sphere gases and sohds. 
Only rarely are exact solutions of these equations possible. Probably the most important technique for generating 
approximate solutions to one-body kinetic equations is the Chapman-Enskog method which, as originally formulated, 
consists of a gradient expansion about a local-equilibrium state The goal in this approach is to construct 
a particular type of solution, called a "normal solution", in which all space and time dependence of the one-body 
distribution occurs implicitly via its dependence on the macroscopic hydrodynamic fields. The latter are, for a 
simple fluid, the density, velocity and temperature fields corresponding to the conserved variables of particle number, 
momentum and energy respectively. (In a multicomponent system, the partial densities are also included.) The 
Chapman-Enskog method proceeds to develop the solution perturbatively in the gradients of the hydrodynamic fields: 
the distribution is developed as a functional of the fields and their gradients and at the same time the equations 
of motion of the fields, the hydrodynamic equations, are also developed. The zeroth-order distribution is the local- 
equilibrium distribution; at first order, this is corrected by terms involving linear gradients of the hydrodynamic 
fields which in turn are governed by the Euler equations (with an explicit prescription for the calculation of the 
pressure from the kinetic theory). At second order, the hydrodynamic fields are governed by the Navier-Stokes 
equations, at third order, by the Burnett equations, etc. The calculations involved in extending the solution to 
each successive higher order are increasingly difficult and since the Navier-Stokes equations are usually considered 
an adequate description of fluid dynamics, results above third order (Burnett order) for the Boltzmann equation 
and above second (Navier-Stokes) order for the Enskog equation are sparse. The extension of the Chapman-Enskog 
method beyond the Navier-Stokes level is, however, not physically irrelevant since only by doing so is it possible to 
understand non-Newtonian viscoelastic effects such as shear thinning and normal stresses which occur even in simple 
fluids under extreme conditions 5], 6]. 

Recently, interest in non-Newtonian effects has increased because of their importance in fluidized granular materials. 
Granular systems are composed of particles - grains - which lose energy when they collide. As such, there is no 
equilibrium state - an undriven homogeneous collection of grains will cool continuously. This has many interesting 
consequences such as the spontaneous formation of clusters in the homogeneous gas and various segregation phenomena 
in mixtures[3|,0,|§|,0|- The coUisional cooling also gives rise to a unique class of nonequilibrium steady states due 
to the fact that the cooling can be balanced by the viscous heating that occurs in inhomogeneous flows. One of the 
most widely studied examples of such a system is a granular fluid undergoing planar Couette flow where the velocity 
field takes the form v (r) = ayx, where a is the shear rate. The common presence of non-Newtonian effects, such 
as normal stresses, in these systems has long been recognized as signalling the need to go beyond the Navier-Stokes 
description jllj. As emphasized by Santos et al,.12], the balance between the velocity gradients, which determine 
the rate of viscous heating, and the cooling, arising from a material property, means that such fluids are inherently 
non-Newtonian in the sense that the sheared state cannot be viewed as a perturbation of the unsheared, homogeneous 
fluid and so the usual Navier-Stokes equations cannot be used to study either the rheology or the stability of the 
sheared granular fluid. One of the goals of the present work is to show that a more general hydrodynamic description 
can be derived for this, and other flow states, which is able to accurately describe such far-from-equilibrium states. 
The formalism developed here is general and not restricted to granular fluids although they do provide the most 
obvious application. Indeed, an application of this form of hydrodynamics has recently been presented by GarzollSj 
who studied the stability of a granular fluid under strong shear. 

The extension of the Chapman-Enskog method to derive the hydrodynamics for fluctuations about an arbitrary 
nonequilibrium state might at first appear trivial but in fact it involves a careful application of the ideas underlying 
the method. To illustrate, let / (r, v, t) be the probability to find a particle at position r with velocity v at time t. 
For a D— dimensional system in equilibrium, this is just the (space and time-independent) Gaussian distribution 

/(r,v,i)=^o(v;n,r,C/)=n(^^-^j exp (- (v - U) V^^bT) (1) 

where n is the number density, kg is Boltzmann's constant, T is the temperature, m is the mass of the particles and 
U is the center-of-mass velocity. The zeroth-order approximation in the Chapman-Enskog method is the localized 
distribution /W (r, v, t) = 0o (v; n (r, t) , T (r, i) , U (r, t)) or, in other words, the local equilibrium distribution. In 
contrast, a homogeneous non-equilibrium steady state might be characterized by some time-independent distribution 

/(r,v,i) = $,3(v;n,T,U) (2) 
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but the zeroth-order approximation in the Chapman-Enskog method will not in general be the localized steady-state 
distribution, /("^ (r,v,t) ^ (v; n (r, t) , T (r, i) , U (r, i)). The reason is that a steady state is the result of a balance 
- in the example given above, it is a balance between viscous heating and collisional cooling. Thus, any change in 
density must be compensated by, say, a change in temperature or the system is no longer in a steady state. This 
therefore gives a relation between density and temperature in the steady state, say n = nss{T), so that one has 
^ss (v;7T,,r, U) = (v; nss(T), r, U). Clearly, it makes no sense to simply "localize" the hydrodynamic variables 
as the starting point of the Chapman-Enskog method since, in a steady state, the hydrodynamic variables are not 
independent. Limited attempts have been made in the past to perform the type of generalization suggested here. 
In particular, Lee and Dufty considered this problem for the specific case of an ordinary fluid under shear with an 
artificial thermostat present so as to make possible a steady state Q,^. However, the issues discussed in this paper 
were circumvented through the use of a very particular type of thermostat so that, while of theoretical interest, that 
calculation cannot serve as a template for the more general problem. 

In Section II, the abstract formulation of the Chapman-Enskog expansion for fluctuations about a non-equilibrium 
state is proposed. It not only requires care in understanding the zeroth order approximation, but a generalization in 
the concept of a normal solution. In Section III, the method is illustrated by application to a simple kinetic theory for 
a sheared granular gas. Explicit expressions are given for the full complement of transport coefficients. One unique 
feature of the hydrodynamics obtained in this case is that several transport coefficients depend linearly on fluctuations 
in the velocity in the y-direction (i.e. in the direction of the velocity gradient). The section concludes with a brief 
summary of the resulting hydrodynamics and of the linearized form of the hydrodynamic equations which leads to 
considerable simplification. The paper ends in Section IV with a summary of the results, a comparison of the present 
results to the results of the standard Chapman-Enskog analysis and a discussion of further applications. 



II. THE CHAPMAN-ENSKOG EXPANSION ABOUT AN ARBITRARY STATE 



A. Kinetic theory 

Consider a single-component fluid composed of particles of mass m in D dimensions. In general, the one-body 
distribution will obey a kinetic equation of the form 

^+v-V^/(r,v,t) = J[r,v,i|/] (3) 

where the collision operator J[r,v,t\f] is a function of position and velocity and a functional of the distribution 
function. No particular details of the form of the collision operator will be important here but all results are formulated 
with the examples of BGK-type relaxation models, the Boltzmann equation and the Enskog equation in mind. The 
first five velocity moments of / define the number density 

n(r,t)= J dv/(r,v,i), (4) 



the flow velocity 



and the kinetic temperature 



n{v,t)^-^j d^^f{v,^,t), (5) 



nr,t) ^ —^^^ [ dvC\r,t)f{r,v,t), (6) 
Dn(r,t)kB J 

where C(r,t) = v — u(r, t) is the peculiar velocity. The macroscopic balance equations for density n, momentum mu, 
and energy i^nksT follow directly from eq. Q by multiplying with 1, mv, and ^mC^ and integrating over v: 

An + nV-u (7) 

DtUi + {mny^V jPij = 
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where Dt = i9t + u • V is the convective derivative. The microscopic expressions for the pressure tensor P = P [/], the 
heat flux q = q [/] depend on the exact form of the coUision operator (see ref. for a general discussion) but as 

indicated, they are in general functionals of the distribution, while the cooling rate ^ is given by 

= Dnir,t)ksTir,t) j ^^^^^t^^ v, t|/]. (8) 



B. Formulation of the gradient expansion 

The goal of the Chapman-Enskog method is to construct a so-called normal solution to the kinetic equation, eq.|(21l. 
In the standard formulation of the method 1], this is defined as a distribution /(r,v,t) for which all of the space 
and time dependence occurs through the hydrodynamic variables, denoted collectively as tp = {n, u, T}, and their 
derivatives so that 

/(r, V, = / (v; i; (r, t) , W (r, t) , VV^ (r, t) ,...)• (9) 

The distribution is therefore a functional of the fields ?/; (r, t) or, equivalently in this case, a function of the fields 
and their gradients to all orders. In the following, this particular type of functional dependence will be denoted more 
compactly with the notation / (v; [V^"-* (r, t)] ) where the index, n, indicates the maximum derivative that is used. 
When all derivatives are possible, as in eq.© the notation /(r, v, i) = / (v; [V'^°°)'0 (r, t)\ ) will be used. The kinetic 
equation, eq.|(21), the balance equations, eqs.Q, and the definitions of the various fluxes and sources then provide a 
closed set of equations from which to determine the distribution. Note that since the fluxes and sources are functionals 
of the distribution, their space and time dependence also occurs implicitly via their dependence on the hydrodynamic 
fields and their derivatives. 

Given such a solution has been found for a particular set of boundary conditions yielding the hydrodynamic state 
-00 (r, t) with distribution /o (v; [V^°°-''0o (r, t)\ ) , the aim is to describe deviations about this state, denoted 54'^ so that 
the total hydrodynamic fields are if) = il)^^ + Sip. In the Chapman-Enskog method, it is assumed that the deviations 
are smooth in the sense that 

5V > ^Vf^V- > ^^VVJV'..., (10) 

where / is the mean free path, so that one can work perturbatively in terms of the gradients of the perturbations to 
the hydrodynamic fields. To develop this perturbation theory systematically, it is convenient to introduce a fictitious 
small parameter, e, and to write the gradient operator as V = V'"' -I- eV*-^'' where the two operators on the right are 
defined by Vo'0 Vi/^o and Vitjj — VStjj. This then generates an expansion of the distribution that looks like 

/(v;[v(°°)V(r,t)]) = /(°)(v;V^°°V(r,i)) (11) 

-fe/W (v;Vi<5V',v(°°V(r,t)) 

+£2/(2) (v; Vi Vi^^, {V,6i;f , V^°^V (r, t)) 
+... 

where /(^^ will be linear in WiSip, Z^^) will be linear in ViVi(5'0 and (WiSip)'^, etc. This notation is meant to be 

taken literally: the quantity Vq°°V (r, t) — {ip (r, t) , VoV' (r, t) , ...} = {^p (r, t) , V-fAo (r, t) , ...} so that at each order in 
perturbation theory, the distribution is a function of the exact field ip (r, t) as well as all gradients of the reference 
field. This involves a departure from the usual formulation of the Chapman-Enskog definition of a normal state. 
In the standard form, the distribution is assumed to be a functional of the exact fields ip (r, t) whereas here it is 
proposed that the distribution is a functional of the exact field tp {r,t) and of the reference state ipo (r, t). Of course, 
it is obvious that in order to study deviations about a reference state within the Chapman-Enskog framework, the 
distribution will have to be a functional of that reference state. Nevertheless, this violates, or generalizes, the usual 
definition of a normal solution since there are now two sources of space and time dependence in the distribution: 
the exact hydrodynamics fields and the reference hydrodynamic state. For deviations from an equilibrium state, this 
point is moot since Vf/'o (r, t) = 0, etc. 

The perturbative expansion of the distribution will generate a similar expansion of the fluxes and sources through 
their functional dependence on the distribution, see e.g. eq.ljSJ), so that one writes 

Pv-P}P+^P}'^ + - (12) 
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and so forth. Since the balance equations Hnk space and time derivatives, it is necessary to introduce a multiscale 
expansion of the time derivatives in both the kinetic equation and the balance equations as 



= + + (13) 

The precise meaning of the symbols dj:^\ d^'^^ is that the balance equations define dj:^^ in terms of the spatial gradients 
of the hydrodynamic fields and these definitions, together with the normal form of the distribution, define the action 
of these symbols on the distribution. Finally, to maintain generality, note that sometimes (specifically in the Enskog 
theory) the collision operator itself is non-local and must be expanded as well in gradients in Sip so that we write 

J[r,v,i|/] = Jo[r,v,i|/]+eJi[r,v,i|/] + ... (14) 

and it is understood that Jo[r, v, t\f] by definition involves no gradients with respect to the perturbations Sip (r, t) but 
will, in general, contain gradients of all orders in the reference fields tpo {r,t). (Note that the existence of a normal 
solution is plausible if the spatial and temporal dependence of the collision operator is also normal which is, in fact, 
generally the case. However, for simplicity, no effort is made here to indicate this explicitly.) A final property of the 
perturbative expansion concerns the relation between the various distributions and the hydrodynamic variables. The 
zeroth order distribution is required to reproduce the exact hydrodynamic variables via 

n{r,t) \ /• / 1 \ , , 
n(r,t)u(r,i) = / V /W(v;V(°°V(r,i))dv (15) 
Dn{v,t)kBT j \ mC^ J ^ ' 

while the higher order terms are orthogonal to the first three velocity moments 

V I /^"' fv;V^°"V(r,i))dv = 0, ?i>0, (16) 
so that the total distribution / ^ + /(^^ + ... satisfies cqs.^-®. 



C. The reference state 



Recall that the goal is to describe deviations from the reference state i/'o (r, t) which corresponds to the distribution 
/o (r, v,t; [V'o]) and in fact the distribution and fields are related by the definitions given in eqs.Q)-®. The reference 
distribution is itself assumed to be normal so that the dependence on r and t occurs implicitly through the fields. In 
terms of the notation used here, the reference distribution satisfies the kinetic equation, eq.0, and the full, nonlinear 
balance equations, eqs.(J7|). Using the definitions given above, these translate into 

(Qf) + V . V(°)) /o (r, v,i; [^o]) = Jo[r, v, t\h\ (17) 

and the fields are solutions to the full, nonlinear balance equations 

a^°^no+u-V(")no + 7ioV(") • uo 
9("^.o.+Uo-V(°).o., + (mno)-^af /^f ' 

where, e.g., P^'^^ is the pressure tensor evaluated in the reference state, and 

a}"Vo = 0, 71 > 0. (19) 

Thus, in the ordering scheme developed here, the reference state is an exact solution to the zeroth order perturbative 
equations. 

For the standard case describing deviations from the equilibrium state, the hydrodynamic fields are constant in 
both space and time and C'-*"^'' = so that the balance equations just reduce to df^^l^^ = 0. The left hand side of the 



= (18) 

= 
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kinetic equation therefore vanishes leaving = Jo[r,v,t\fo] which is indeed satisfied by the equilibrium distribution. 
For a granular fluid, (^^^^ ^ and the simplest solution that can be constructed consists of spatially homogeneous, 
but time dependent fields giving 



afVo(r,v,i;[^o]) = Jo[r,v,i|/o] (20) 



and 



di°^nn = (21) 







so that the distribution depends on time through its dependence on the temperature. The balance equations, together 
with the assumption of normality, serve to define the meaning of the left hand side of eq. (|20|l giving 

-C^"°)To^/o(r,v,i;[V'o]) = Jo[r, v, t|/o]. (22) 



Typically, this is solved by assuming a scaling solution of the form /o (r, v,t; ["^o]) = $ ( v 



feBT(t) 



D. The zeroth order Chapman-Enskog solution 

As emphasized above, the Chapman-Enskog method is an expansions in gradients of the deviations of the hy- 
drodynamic fields from the reference state. Using the ordering developed above, the zeroth order kinetic equation 
is 

af/'"' (r, v; (r, t) , [V-o]) + v • V^")/^ (r, v; (r, t) , [^o]) = Jo[r, v, t\fo]. (23) 
and the zeroth order balance equations are 

af^n + u-Vno+nV-Uo = (24) 
5f ^M,+u-Vuo.^ + (mn)Tiv(°)/^^°^ = 

a^°)r + u.vro + -^(v(").qW + /^f 9,^.0,.) = -C^'^t. 

Making use of the balance equations satisfied by the reference fields, (|18fl , this can be written in terms of the deviations 

as 

di°hn + Su-V no + 6n\7 -uo = (25) 

Dukb V 'J / DuoKb ^ i-j ■• ; 

Since the zeroth-order distribution is a junction of but a functional of the reference fields, the time derivative in 
ea. (P5)l is evaluated using 

and these equations must be solved subject to the additional boundary condition 

hm (r, V, t; (r, t) , [^o]) = /o (r, v, t; [^o]) • (27) 

There are several important points to be made here. First, it must be emphasized that the reference fields ipo (r, t) and 
the deviations Sip (r, t) are playing different roles in these equations. The former are fixed and assumed known whereas 
the latter are independent variables. The result of a solution of these equations will be the zeroth order distribution 
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as a function of the variables 6ip. For any given physical problem, the deviations will be determined by solving the 
balance equations, eas. H25|l . subject to appropriate boundary conditions and only then is the distribution completely 
specified. Second, nothing is said here about the solution of eqs.(j2SJ-(I2SJ which, in general, constitute a comphcated 
functional equation in terms of the reference state variables ipo^a (r, t) . The only obvious exceptions, and perhaps 

the only practical cases, are when the reference state is either time- independent, so that 9j°''V'o.q = 0, or spatially 
homogeneous so that f^^Hs a function, and not a functional, of the reference fields. The equilibrium state is both, the 
homogeneous cooling state is a spatially homogeneous state and time-independent flow states such as uniform shear 
flow or Pouseille flow with thermalizing walls are important examples of time- independent, spatially inhomogeneous 
states. Third, since eas. (|^ - (l^ are the lowest order equations in a gradient expansion, they are to be solved for 
arbitrarily large deviations of the fields. Sip. There is no sense in which the deviations should be considered to be 
small. The fourth observation, and perhaps the most important, is that there is no conceptual connection between the 

zeroth order distribution /(°) ^v; dtjj (r, t) , Vq°°V-'o (r, t)^ and the reference distribution /o (v; Vo°°Vo (r, t) j except 
for the limit given in ea. (|27|l . In particular, it will almost always be the case that 

(v; (r, t) , V(°°) Vo (r, t)) ^ /„ (v; v(°°) (^o (r, t) + (r, t))) . (28) 

A rare exception for which this inequality is reversed is when the reference state is the equilibrium state. In that case, 
the density, temperature and velocity fields are uniform and the reference distribution is just a Gaussian 

/o(r,v;V^°°Vo) = 0o (v; no. To, Uo) (29) 
and the solution to the zeroth order equations is the local equilibrium distribution 

/(°) (v; JtA (r, t) , v[,°°Vo (r, t)) = (v; n + ,5n (r, t) , T + ST (r, t) , U + 5U (r, t)) - /o (v; V^°") (^o (r, t) + S^j (r, i))) . 

(30) 

For steady states, as will be illustrated in the next Section, it is not the case that /^°^ is obtained from the steady-state 
distribution via a "localization" along the lines of that shown in ea. (|30|l . On the other hand, eas. H23|) - (|24|l are the 
same whether they are solved for the general field Sip (r, t) or for the spatially homogeneous field Sip (t) with the 
subsequent localization Sip (t) — > S'4>{r,t). Furthermore, these equations are identical to those one would solve in 
order to obtain an exact normal solution to the full kinetic equation, ea. (|17|l and balance equations, ea. H18() . for the 
fields ipo (r, t) + Sip (t). In other words, the zeroth-order Chapman-Enskog distribution is the localization of the exact 
distribution for homogeneous deviations from the reference state. Again, only in the case of the equilibrium reference 
state is it true that this corresponds to the localization of the reference state itself. 



E. First order Chapman-Enskog 

In the following, the equations for the first-order terms will also be needed. Collecting terms in ea. p7|l . the first 
order distribution function is found to satisfy 

S'^V^'^v; Sip (r, t) , [iPo]) + V • VW/(i)(v; SiJj (r, t) , [^o]) (31) 

= Jo[r, V, t\h] + Ji [r, V, t\fo] ~ [d^'^ f^°^ (v; SiP (r, t) , [V-o]) + v • V(i)/(°)(v; SiP (r, t) , [V-o])) 

and the first-order balance equations become 

d^t^^Sn + u-VSn + nV -Su (32) 
dPSu. + u-WSui + imnj-^wf^P^f + {mny^wfpj;p = 



III. APPLICATION TO UNIFORM SHEAR FLOW OF GRANULAR FLUIDS 



Uniform shear fiow (USF) is a macroscopic state that is characterized by a constant density, a uniform temperature 
and a simple shear with the local velocity field given by 

Ui = a^jrj^ — (xS'ixSjy^ (33) 
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where a is the constant shear rate. If one assumes that the pressure tensor, heat flux vector and coohng rate are also 
spatially uniform, the reference-state balance equations, eas. (|18ll . become 



d^"^no - (34) 
df°^uo^i + aua^ydix = 

The question of whether or not these assumptions of spatial homogeneity are true depends on the detailed form 
of the collision operator: in ref.|l7| it is shown that only for the linear velocity profile, ea. (|33|) . this assumption is 
easily verified for the Enskog kinetic theory (and hence for simpler approximations to it such as the Boltzmann 
and BGK theories). This linear velocity profile is generated by Lee-Edwards boundary conditions which are 
simply periodic boundary conditions in the local Lagrangian frame. For elastic gases, ^C^^) = and the temperature 
grows in time due to viscous heating and so a steady state is not possible unless an external (artificial) thermostat 
is introduced|l5l|. However, for inelastic gases, the temperature changes in time due to the competition between 
two (opposite) mechanisms: on the one hand, viscous (shear) heating and, on the other hand, energy dissipation in 
collisions. A steady state occurs when both mechanisms cancel each other at which point the balance equation for 
temperature becomes 

-^aPi°°) = -C(°°)To. (35) 

Note that both the pressure tensor and the cooling rate are in general functions of the two control parameters, the 
shear rate and the coefficient of restitution, and the hydrodynamic variables, the density and the temperature, so 
that this relation fixes any one of these in terms of the the other three: for example, it could be viewed as giving the 
steady-state temperature as a function of the other variables. 

At a microscopic level, the one-body distribution for USF will clearly be inhomogeneous since the eq.© and ea. H33fl 
imply that the steady-state distribution must give 

ayx=— / fivv/o(r,v). (36) 
no J 

However, it can be shown, at least up to the Enskog theory that for the Lees-Edwards boundary conditions, the 
state of USF possesses a modified translational invariance whereby the steady state distribution, when expressed in 
terms of the local rest- frame velocities Vi = Vi — CLijTj does not have any explicit dependence on position. In terms of 
these variables, and assuming a steady state, the kinetic equation becomes 

-a^,^/(V) = J[V|/,/] . (37) 

The solution of this equation has been considered in some detail for the BGK- type models 0|, ^^,|23|, the 
Boltzmann equation^, and the Enskog equation _21j , .22] , 17.J . 

A. The model kinetic theory 

Here, for simplicity, attention will be restricted to a particularly simple kinetic theory which nevertheless gives 
realistic results that can be compared to experiment. The kinetic theory used is the kinetic model of Brey, Dufty and 
Santos 20], which is a relaxation type model where the operator J[/, /] is approximated as 

J[L /] ^ («) ^ W if - M + \c («) ''W-^- (c/) . (38) 

The right hand side involves the peculiar velocity C = v — u = V — (5u and the local equilibrium distribution, 
eq.(^. The parameters in this relaxation approximation are taken so as to give agreement with the results from the 
Boltzmann theory of the homogeneous cooling state as discussed in ref.[23. Defining the collision rate for clastic hard 
spheres in the Boltzmann approximation as 



-W= (i. + 2)r(i./2) - ^''^ 



8 



the correction for the effect of the inelasticity is chosen to reproduce the Navier-Stokes shear viscosity coefficient of 
an inelastic gas of hard spheres in the Boltzmann approximation . |24l| giving 

'^*i»)^^{'^ + '^)iiD-l)a + D + l). (40) 

The second term in eq.(j2Hl) accounts for the coUisional cooling and the coefficient is chosen so as to give the same 
cooling rate for the homogeneous cooling state as the Boltzmann kinetic theory [ill, nil, 

C* («) = ^ (1 - ■ (41) 

In this case, the expressions for the pressure tensor, heat-flux vector and cooling rate take particularly simple forms 
typical of the Boltzmann description^ 

P^J = mj dCaQ/(r,C,i), (42) 
= J dCQC^f{r,C,t), 
while the cooling rate can be calculated directly from eas. H38f) and JHJl with the result ({''P) — v (ijj) (* (a)- 



B. The steady-state 

Before proceeding with the Chapman-Enskog solution of the kinetic equation, it is useful to describe the steady 
state for which the distribution satisfies eq.||2ZJl which now becomes 

-«n^/(V) = (a) V 0,0) if - .^o) + ^C* («) ^ (V-o) ^ • (V/) . (43) 
The balance equations reduce to 

2aP^'y = -C* (a) V (^) DnoksTo. (44) 
An equation for the pressure tensor is obtained by multiplying eg . (|43|l through by raViVj and integrating giving 

aP-y5,^ + aP;^5,^ = -v* (a) v (V-o) (i^7 - nofesToJ^ ) - C* («) v (V'o) Pt; ■ 
This set of algebraic equations is easily solved giving the only non-zero components of the pressure tensor as 

P- - -*H+^-ff(") „„fc,To (45) 
V* (a) -I- C* (a) 

pas p 

j.*(a)+C(a) 

where a*^ — agsjv (V'o) satisfies the steady-state condition, ea. (|44|l 

5£ ^ Q, . (46) 

For fixed control parameters, a and a, this is a relation constraining the state variables no and Tg. The steady-state 
distribution can be given explicitly, see e.g. |25j |. 

C. Zeroth order Chapman-Enskog 

Since the only spatially varying reference field is the velocity and since it is linear in the spatial coordinate, the 
zeroth-order kinetic equation, eq.JSSJ becomes 

+ V . (v(")uo.) ^/(") = -u W (/(") - 0o) + IC {a) K^) I; • (c/(°)) • (47) 
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or, writing this in terms of the pecuhar velocity, 



9f /("^ + - avy^j^"^ ^ -V (^) (/ - 0o) + \C (a) W |^ • (c/(«)) . (48) 

Here, the second term on the left accounts for any explicit dependence of the distribution on the coordinate y, aside 
from the implicit dependence coming from C. Since it is a zero-order derivative, it does not act on the deviations Si/j. 
In terms of the peculiar velocity, this becomes 

5f + {Cy + suy) dif'^ - -Cy^/"' - ""^""yw:/'' = - + V* " A ■ • 

The first term on the left is evaluated using ea. (|26|l and the zcroth order balance equations 

= (50) 
df^u^ + aSuySix = 

and the assumption of normality 

arf , ^ {a?>^n) (^/<«.) + {af'^T) (^/<».) + {a<%u,) (^/<». 

to give 

= (a) V (,,!.) (/<»' - *„) + {o.) (VI ^ ■ (c/<»') , 

where the temperature derivative is understood to be evaluated at constant density. Here, the second term on the 
left in eq. (|49|l has been dropped as neither eq. H49|l nor the balance equations contain explicit reference to the velocity 
field uo, and so no explicit dependence on the coordinate y , thus justifying the assumption that such dependence 
does not occur in /(o). One can also assume that /^°) depends on Sui only through the peculiar velocity, since in that 
case the term proportional to ^Uj^ vanishes as well and there is no other explicit dependence on 6uy. 

Equation (|51|l is closed once the pressure tensor is specified. Since the primary goal here is to develop the transport 
equations for deviations from the reference state, attention will be focused on the determination of the pressure tensor 
and the heat flux vector. It is a feature of the simple kinetic model used here that these can be calculated without 
determining the explicit form of the distribution. 



1. The zeroth-order pressure tensor 



An equation for the pressure tensor can be obtained by multiplying this equation through by mCiCj and integrating 
over velocities. Using the definition given in ea. H42|) . 

(^-C («) ^ W T ^^aPi^'^ ^P,f + a<5„4°) + aS^^P^^ = -v* {a) v (V^) (/^f - k^nksT) - C (a) v (^) p\f , 

(52) 

and of course there is the constraint that by definition Tr(P) — DnksT . It is interesting to observe that eqs.,|SD| - 

(|52|l are identical with their steady-state counterparts when the steady-state condition, C}-'^^T = j^^^ aPxy , is fulfilled. 
However, here the solution of these equations is needed for arbitrary values of i5T, 8n and (5u. Another point of interest 
is that these equations are local in the deviations so that they are exactly the same equations as describe spatially 
homogeneous deviations from the reference state. As mentioned above, this is the meaning of the zeroth-order solution 
to the Chapman-Enskog expansion: it is the exact solution to the problem of uniform deviations from the reference 
state. It is this exact solution which is "localized" to give the zeroth-order Chapman-Enskog approximation and not 
the reference distribution, /o, except in the rare cases, such as equilibrium, when they coincide. 
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To complete the specification of the distribution, eqs. (|51|l and H52() must be supplemented by boundary conditions. 
The relevant dimensionless quantity characterizing the strength of the nonequilibrium state is the dimensionless shear 
rate defined as 



iD + 2)r{D/2) ma^ 

It is clear that for a uniform system, the dimensionless shear rate becomes smaller as the temperature rises so that 
we expect that in the limit of infinite temperature, the system will behave as an inelastic gas without any shear - i.e., 
in the homogeneous cooling state, giving the boundary condition 



1 

T— foo nksT 



lim = (54) 



and in this limit, the distribution must go to the homogeneous cooling state distribution. These boundary conditions 
can be implemented equivalently by rewriting eas. HftUI) in terms of the inverse temperature, or more physically the 
variable a* , and the dimensionless pressure tensor P>- ' — _^^^^ py ' giving 

and writing /(o) (C; V) = n ""^^ 9 {/j 



D/2 



—V* (a) (g — exp [—mC^/kBTy) , 



(56) 



with boundary condition lim^'^o Pij'^ — ^ij ^^^d lim^.^g g = exp (^—mC'^/kBT^- For practical calculations, it is more 
convenient to introduce a fictitious time variable, s, and to express these equations as 

^ = \a*Q*{a) + ^a*'Pi*J (57) 

where the boundary condition is then P^^j"^ (s = 0) = 5ij^ and a* (s = 0) = 0. The distribution then satisfies 

^5 = -^a*P^*y>a—g - a*P^;^g + a*Cy—g - v* (a) {g - exp {-mC^/kBT)) (58) 

with lims^o 9 — exp (^—mC^ /kBPy These are to be solved simultaneously to give Pj-*^ (s) , a* (s) and f^^^ (s) from 

which the desired curves Pj^*^ (a*) and /^°^ (a*) are obtained. 

Physically, if the gas starts at a very high temperature, it would be expected to cool until it reached the steady 
state. It is easy to see that the right hand sides of eas. (|57|l do in fact vanish in the steady state so that the steady 
state represents a critical point of this system of differential equationsl^^. In order to fully specify the curve Pij (T) 
and the distribution /(°) it is necessary to integrate as well from a temperature below the steady state temperature. 
Clearly, in the case of zero temperature, one expects that the pressure tensor goes to zero since this corresponds to the 
physical situation in which the atoms stream at exactly the velocities predicted by their positions and the macroscopic 
flow field. (Note that if the atoms have finite size, this could still lead to collisions. However, the BGK kinetic theory 
used here is properly understood as an approximation to the Boltzmann theory appropriate for a low density gas in 
which the finite size of the grains is of no importance.) Thus, the expectation is that the zero-temperature limit will 
give 

lim Pjf^ = 0. (59) 

Then, in terms of a fictitious time parameters, one has 

rIT 9 

— ^ -C{a)u{^)T--aTPi*^ (60) 
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and for the distribution 



A final point is that the solution of these equations requires more than the boundary condition P^^^ (s — 0) — since 



evaluation of the right hand side of ea. l|()U|) requires a statement about Pj;*'' (s = 0) as well. A straight-forward series 



solution of ea. (|52|l in the vicinity of T = gives P*y ^ a* -^^^ and P*^ ^ a* ^/■^so that the correct boundary condition 

here is pj^*^ (s = 0) — 0. The solution of these equations can then be performed as discussed in ref. |l,2j with the 
boundary conditions given here. 

It will also prove useful below to know the behavior of the pressure tensor near the steady-state. This is obtained 
by making a series solution to ea. (|55|l in the variable (a* — a*^) where a*^ is the reduced shear in the steady-state. 
Details are given in Appendix 1X1 and the result is that 



= 1 + 4*, )+•••)' (62) 



with the coefficients 

A( a) + C («) 
C («) 



\A{a) + iy* (a) + [a) J 
( ^ (A (a) + j^v* [a) + \C ia)) jy* («) + C («)) 

(A (a) + V* (a) + \C (a)) {v* (a) + DC (a)) ' 

where A (a) is the real root of 

{y* (a) C* (a)) ^^^{^v*'^ (a) -h 14:/* (a) C* (a) + 7C*^ (a)) A-^C* (a) (2^*^ (q,) _ ^* (q,) _ 2^*2 ^^^^ ^ q_ 

(64) 



2. Higher order moments: the zeroth-order heat flux vector 

Determination of the heat flux vector requires consideration of the full tensor of third order moments. Since fourth 
order moments will also be needed later, it is easiest to consider the equations for the general A^-th order moment, 
defined as 

Kl.N (r,i) = mjd^, Q,...C.„/(") (r, C,i) . (65) 

To simplify the equations, a more compact notation will be used for the indices whereby a collection of numbered 
indices, such as will be written more compactly as In so that capital letters denote collections of indices and 

the subscript on the capital indicates the number of indices in the collection. Some examples of this are 

<^ = (66) 

-12 Il«2 

I2V 1112V 



In terms of the general moments, the heat flux vector is 



) (r,i) = i ^ M^g (r,t) = \mI]] (r,t) , (67) 



where the second equality introduces the Einstein summation convention whereby repeated indices are summed. The 

0(0) ^ . AO] 



pressure tensor is just the second moment P^'^^ = ^^ij^ ■ The local equilibrium moments are easily shown to be zero 



for odd A^ while the result for even A^ is 

Me 
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where the operator Vijk... indicates the sum over distinct permutations of the indices ijk... and has no effect on any 
other indices, (e.g., Vi^Si^i^Si^i^ = 6i-^i^di^i^ + Si-^i^Si^i^ + Si^i^Si^i^). An equation for the general A^-th order moment 
can be obtained from ea. l47|l with the result 

(-C* («) - T-^mI'J + (^^* (a) + («)) mI°J + a*ViJ..,Mfl^y = i.* (a) uf^^^. (69) 

Writing Mfj = mn {^) ' Af/^gives 

- (^C (a) + |a*Pi;)) r^M;„ + (^^* (a) - + aViJ.^,Ml_^y = (a) Nlf^^*^ (70) 

Notice that the moments are completely decoupled order by order in N . Since the source on the right vanishes 
for odd N it is natural to assume that Mj^ ~ for odd N . This is certainly true for temperatures above the 
steady state temperature since the appropriate boundary condition in this case, based on the discussion above, is that 
limr^oo = ^f^*^ ~ 0. In the opposite limit, T ^ 0, as mentioned above, one has that P*y ~ a*~^/^ ~ 7^1/6 
and there are two cases to consider depending on whether or not the third term on the left contributes. If it does, 
i.e. if one or more indices is equal to x, then a series solution near T = gives M;^ ~ a*-i - while if no index 
is equal to x then Mj^ ~ a*~'^/^ ~ 7^i/3 giving in both cases the boundary condition limT^o = 0- In particular, 
this shows that the odd moments vanish for all temperatures. From this, it immediately follows that 

) {v,t) = 0. (71) 

D. First-order Chapman-Enskog: General formalism 

The equation for the first-order distribution, ea. H31|) . becomes 

d['^f('^ + avy^/'^ = -V W /(^) + IC {a) H^)^- (c/(^)) - (af + V • Vi/(")) , (72) 

and the operator d^^^ is defined via the corresponding balance equations which are now 

dl^^Sn + u-V6n + nV ■ 5u ^ (73) 
5t'^'(5w, + u- V(5ii, + (mn)-iaf^p/°^ + (mn)-i9(°)i^(J^ = 



Writing the kinetic equation in the form 

d 



dl"^f'^ + ^^Jvf^'^ + ^* (") ^ W f^'^ - IC («) ^W-^- (C/^'^) (74) 

equations for the A^-th moment can be obtained by multiplying through by Ci-^...Cij^ and integrating over velocity. 
The first two terms on the left contribute 
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where the last hne follows from using tli6 z6rot}i order balance equcition '^Suij^ — ~^^i^x 

5uy . The evaluation of 

the right hand side is straightforward with the only difficult term being 



(76) 



and from ea. (|7U|l it is clear that Mj-J is independent of Suj so that the first term on the right vanishes. Thus 

di°h-4'J + (m«^ + <5u,M«) +aViJ..,Mll\^^ + (^i.* (a) + («)) ^ W m'^'J (77) 

- (din,) m["J - Vj^ Mfl^, - dlM^l] 



Superficially, it appears that the right hand side depends explicitly on the reference field, since ui — Uqj + Sui, which 
would in turn generate an explicit dependence of the moments on the y-coordinate. However, when the balance 
equations are used to eliminate d^^^ this becomes 



9f Mj^) 

t -IN 



+ «^ {m[1\ + SuyM['J) + aViJ^^^M^l^^ + (v* (a) + («)) W M<i'J (78) 



dT 



M 



(0) 



Then, assuming that the first-order moments are independent of the reference field, Uo, gives 



I) { 2a d 



DnkB dT I 



d 



(0) 



(0) 



Dnk 



B 



dT 



+ 



In-1 Q§j^ InI 



--Pi, 



1 



_^p(0)\,^(0) 



{d^^^ST) 



(79) 



which is consistent since no factors of Uq appear and since the zeroth order moments are known to be independent of 
the reference velocity field. 

The moment equations are linear in gradients in the deviation fields, so generalized transport coefficients via the 
definition 



M 



(1) 



-A 



In ab ' 



d&tpb 
dra 



dSn dST dSu^ 

VlNab- 



On 



(80) 



where the transport coefficients for different values of N have the same name but can always be distinguished by the 
number of indices they carry. The zeroth-order time derivative is evaluated using 



dra 



d\lf,ab\ dSlpb 



d\">T 



dT 
dT dra 



+ {di"'Su,) 



(81) 



dSuj J dra 



dSip, 



d 



{di'^S^b) 



dra dSipc 



(82) 



where the third line follows from (a) the fact that the transport coefficients will have no explicit dependence on the 
velocity field, as may be verified from the structure of ea H79|) and (b) the fact that the gradient here is a first order 
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gradient Vi so that it only contributes via gradients of the deviations of the fields thus giving the last term on the 
right. Since the fields are independent variables, the coefficients of the various terms must vanish independently. 
For the coefficients of the velocity gradients, this gives 

/ f \ \ d ^ ('9t''''<5'"c ) / j\f \ f 2a d ( \\ 

[dt'Tj g^ViNub + ViNac + aVi^6^i^r]i^_,yab + (a) + yC* (a) j i' (V') Vinab - ( ) 



^ 'r(o)_nAo)\ 2 (0) 9 (0) (o) 



The vanishing of the coefBcients of the density gradients gives 



mn^" U<5n dSn 



while the vanishing of the coefficient of the temperature gradient gives 



Notice that for even moments, the source terms for the density and temperature transport coefficients all vanish (as 
they involve odd zeroth-order moments) and it is easy to verify that the boundary conditions are consistent with 
Hiiva = K/jva = and only the velocity gradients contribute. For odd values of N, the opposite is true and Tj/^aft = 
while the others are in general nonzero. 

E. Navier-Stokes transport 

1. The first order pressure tensor 

Specializing to the case N = 2 gives the transport coefficients appearing in the pressure tensor 



(1) dSua 



'In =-^lNab-^ (86) 



where the generalized viscosity satisfies 

[df^T^ -^Vijab - arjijaxSby + aS^iTjjyab + aS^jTjiyab + {i^* (a) + C* (a)) W Vijab - (^'^^'^^iJ^^ Vxyl§>7) 

2. First order third moments and the heat flux vector 

For the third moments, the contribution of density gradients to the heat flux is well-known in the theory of granular 
fluids and the transport coefficient is here the solution of 

^Myte + ^ Kijka + (ly* {a) + (a)) W l^ijka (88) 

mn\dn mn\dn '^^ mn\dn dn 
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and the generalized thermal conductivity is determined from 



(af ^t) —n,,,a + ^ n,,ka + [ly* {a) + -(* (a) j (^) K,,ka (89) 

1 ^ ^ p(o)^ p(o) _ J_ Ap(o)^ p(o) - J- — + AAf(°' 



mn '^'^ y mn V °' 7 ''i mn \dT ) dT 

Note that both of these require knowledge of the zeroth-order fourth velocity moment M^^jl^ . The heat flux vector is 

(1) _ dSn _ d6T 

Ora OTa 

where 

M^a = Myja (91) 

F. The second-order transport equations 

In this Section, the results obtained so far are put together so as to give the Navier-Stokes equations for deviations 
from the steady state. The Navier-Stokes equations result from the sum of the balance equations. To first order, this 
takes the form 

dtn + u- W6n + nW ■ Su = (92) 

where dt = dj:^^ + dj:^\ By analogy with the analysis of an equilibrium system, these will be termed the Euler 
approximation. Summing to second order to get the Navier-Stokes approximation gives 

dtn + u- V5n + nV • (5u = (93) 
dtu, + u • V6u, -f a5,Juy + [mn)-^ df^ P^f + {mn)-^ d\^^ P^l^ + {mn)-^ df^ P^f = 

dtT + vi- VST + (V(i) • q(i) + V(") • q(2) + P^pV,6u, + Pg^VjSu, + aP^^J + aPg^ + aP^f^ = -(* (a) v {i,) T. 

where now dt = df^ + di^^ + df^ but this expression is problematic. Based on the results so far, it seems reasonable to 
expect that dj^^P^p = V^"^ • q*-^-* — 0. However, to consistently write the equations to third order requires knowledge 

(2) 

of Pxy which is not available without extending the solution of the kinetic equation to third order. The reason this 
problem arises here, and not in the analysis about equilibrium, is that the shear rate, a, arises from a gradient of 
the reference field. In the usual analysis, such a term would be first order and aPxy — {diUoj) P^j would be of third 
order and therefore neglected here. This is unfortunate and shows that this method of analysis does not completely 
supplant the need to go beyond the second-order solution in order to study shear flow. However, this problem is not 
unique. In fact, in calculations of the transport coefficients for the homogeneous cooling state of a granular gas, a 
similar problem occurs in the calculation of the cooling rate: the true Navier-Stokes expression requires going to third 
order in the solution of the kinetic eauationp^.p^. (This is because the source does not appear under a gradient, 
as can be seen in the equations above.) Thus, it is suggested here that the same type of approximation be accepted 

(2) 

here, namely that the term aPxy is neglected, so that the total pressure tensor and heat flux vectors are 

P^, = pS^+P^ (94) 

(0) , (1) 
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and the transport equations can be written as 

dtn + V ■ (nu) = (95) 
dtUi + u ■ Vui + {mn)~^djPij ~ 



Unkr 



which is the expected form of the balance equations. The total fluxes are given terms of the generalized transport 
coefficients 



= P^' - (96) 



j(0) ^ d5ua 
dn 

dSn dST 



G. Linearized second-order transport 

Some simplification occurs if attention is restricted to the linearized form of these equations. This is because, as 
noted in the previous Section, several transport coefficients are proportional to 6uy and consequently do not contribute 
when the transport coefhcients are linearized. Taking this into account, the total fluxes are 

p.. = p(^^)+(^] Sn+(^] ST-if^,^ (97) 



dSn ^B5T 



dra '^"^ dra ' 



where the superscript on the transport coefficients, and subscript on the derivatives, indicates that they are evaluated 
to zeroth order in the deviations, 



\ dSn j Stp^Q dSn ' 

i.e. in the steady state. The defining expressions for the transport coefficients simplify since the factor df^T is at least 
of first order in the deviations from the steady state (since it vanishes in the steady state) so that the temperature 
derivative can be neglected thus transforming the differential equations into coupled algebraic equations. Also, all 
remaining quantities are evaluated for = 0, i.e. in the steady state. Thus the viscosity becomes 

-a*ssVtjax5by + aL^^r/l^ab + <s5xjVtyab + («) + C («) ) Vt^ab - ijf^S^ ^ '^-V-b ^^^^ 

where a*^ was defined in ea. H4()l) . The generalized heat conductivities will be given by the simplified equations 



and 



^ ^ka + («) + \C («)) ^ka + «L^..fc4feAC^,a (100) 

'^~'(^o) p ( d (o)\ piss) ,-1 (^i^f d (0) ' 
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In these equations, the hydrodynamic variables 4'Q must satisiy the steady state balance condition, ea. 1)46(1 . The 
various quantities in these equations are known from the analysis of the zeroth order moments. For example, from 
ea. (|H^ . one has that 

-^T,-'P,^fA*^ (a) (101) 



where A*j (a) was given in eg . 1)63(1 and here, there is no summation over repeated indices. The derivatives of higher 
order moments in the steady state can easily be given using the results in Appendix ^The linearized transport 
equations are 

dtSn + ay—5n + tiqV • 5u = (102) 
ox 

„ ™ 9 2 / „„ d'^Sn „„ d'^6T iss)dSui „„ dSua\ 

^DnlkB p 1^ dSn j J ^ Dnoks [ dST J J 

- -^C(a)^(V'o)^T-C(«)^(V^o)ro-. 
2 no 

where the fact that v (ip) ^ nT^/^ has been used. These equations have recently been used by Garzo to study the 
stability of the granular fluid under uniform shear fiow[l3j. 

IV. CONCLUSIONS 

In this paper, the extension of the Chapman-Enskog method to arbitrary reference states has been presented. One of 
the key ideas is the separation of the gradient operator into "zeroth" and "first" order operators that help to organize 
the expansion. It is also important that the zeroth order distribution be recognized as corresponding to the exact 
distribution for arbitrary uniform deviations of all hydrodynamic fields from the reference state. This distribution 
does not in general have anything to do with the distribution in the reference state, except in the very special case 
that the reference state itself is spatially uniform. 

The method was illustrated by application to the paradigmatic non-uniform system of a fluid undergoing uniform 
shear flow. In particular, the fluid was chosen to be a granular fluid which therefore admits of a steady state. The 
analysis was based on a particularly simple kinetic theory in order to allow for illustration of the general concepts 
without the technical complications involved in, e.g., using the Boltzmann equation. Nevertheless, it should be 
emphasized that the difference between the present calculation and that using the Boltzmann equation would be no 
greater than in the case of an equilibrium fluid. The main difference is that with the simplified kinetic theory, it is 
possible to obtain closed equations for the velocity moments without having to explicitly solve for the distribution. 
When solving the Boltzmann equation, the moment equations are not closed and it is necessary to resort to expansions 
in orthogonal polynomials. In that case, the calculation is usually organized somewhat differently: attention is focussed 
on solving directly for the distribution but this is only a technical point. (In fact. Chapman originally developed his 
version of the Chapman-Enskog method using Maxwell's moment equations while Enskog based his on the Boltzmann 
equationj^. The methods are of course equivalent.) 

It is interesting to compare the hydrodynamic equations derived here to the "standard" equations for fluctuations 
about a uniform granular fluid. As might be expected, the hydrodynamic equations describing fluctuations about the 



i^-' (V-o) 



(0)' 



dT 



dP, 



(0)' 



\ dn 



dT 
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state of uniform shear flow are more complex in some ways than are the usual Navier-Stokes equations for a granular 
fluid, but the similarities with the simpler case are perhaps more surprising. The complexity arises from the fact that 
the transport coefficients do not have the simple spatial symmetries present in the homogeneous fluid where, e.g., there 
is a single thermal conductivity rather than the vector quantity that occurs here. However, just as in homogeneous 
system, the heat flux vector still only couples to density and temperature gradients and the pressure tensor to velocity 
gradients so that the hydrodynamics equations, ea. (|95(l . have the same structure as the Navier-Stokes equations for 
the homogeneous system. 

An additional complication in the general analysis presented here is that the zeroth-order pressure tensor and the 
transport coefficients are obtained as the solution to partial differential equations in the temperature rather than as 
simple algebraic functions. This requires that appropriate boundary conditions be supplied which will, in general, 
depend on the particular problem being solved. Here, in the high-temperature limit, the non-equilibrium effects are of 
no importance and the appropriate boundary condition on all quantities is that they approach their equilibrium values. 
Boundary conditions must also be given at low temperature as the two domains are separated by the steady-state 
which represents a critical point. At low temperatures, there are no collisions and no deviations from the macroscopic 
state so that all velocity moments go to zero thus giving the necessary boundary conditions. A particularly simple 
case occurs when the hydrodynamic equations are linearized about the reference state as would be appropriate for a 
linear stability analysis. Then, the transport properties are obtained as the solution to simple algebraic equations. 

A particular simplifying feature of uniform shear flow is that the flow field has a constant first gradient and, as a 
result, the moments do not explicitly depend on the flow field. This will not be true for more complex, nonlinear flow 
fields. However, the application of the methods discussed in Section H should make possible an analysis similar to 
that given here for the simple case of uniform shear flow. 
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APPENDIX A: STEADY-STATE TEMPERATURE LIMIT 



Recall that in the steady state 



p*{ss) ^ y* {a) + hxDC (g) 



p*(ss) 



aiy (a) 



{v* {a)+C* {a)f' 



and the explicit form of the steady-state condition, ea. lfl^ giving the value of the reduced shear in the steady state, 
aL> is 

afsV* («) D 
{ly* (a) + C {a)y 2 

Assume that the stresses are analytic in a* so that near the singularity 

/f)=i>;(-) + ^,, («*-«:.) + ... 

They satisfy ea. (|53|l . 

Qc* (a) + a^^? = - - - ^* («) (^?^ - (A3) 

so 

{^Pir^ + ^<sAxy^ a:,A, = ^Pi^'^P^r' + i<sAxyPt'' + l<sP!^'^ ^ (A4) 
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or 



^atlA^yAj + at, (^i.* (a) + («)) A,, - la^A^yP^^'^ + afj.^A^y + al%,A,y = u* (a) (p^;^*^ - 5^,) (A5) 



In component form this is 



i^a:lA,yAyy + ( (a) + (")) Ayy - ^af,A,yP<^^^ = u* (a) - l) 



1 

1) 



Substituting 



*2 



A^y = DC/a 
{l-Si^)Aii = (l-5i^)DB/al, 



(A6) 



(A7) 



gives 



BC + (a) + ic* («)) B - ^CP<^^) = ^u* (a) (p;(-) - l) 



and the steady-state condition makes this 



' + («) + ^C* («)) C - y:,CP<^^^ + allB = ^u* (a) a^P^f 



(ss) 



D 



BC + (^u* (a) + (a)) B - ^CP<^^^ = ^u* (a) (p<;^) - l) 



The solution is 



C'+[y* {oc) + lc {a)]C + af,B 



Q — p*(ss) 



1 2C - C* (a) 



(A8) 



(A9) 



DC + v*{(x) + \C (a) 

with C being the real root of 

4C^+8 {v* (a) + C* (a)) C^+{Aiy*^ (a) + 14i/* (a) C* (a) + 7C*^ (a)) C+C* (a) (21/*^ (a) - u* (a) (* (a) - 2C*^ (a)) = 
Finally, it is useful to note that the full pressure has the expansion 



Pij = nkBTP*j 

= n ( ksTss + X 



( dkBT \ 
V da* A 



-)( 



(AlO) 
(All) 



This gives 



xy 



P-+n{kBT.,A,+P<-^ (^)^ 
+ xnkBT^s (A^ - 2a*-'P:l''^) + 



P^l + nksTss (^C - 2a*-'P<^A {a* - + ... 



(A12) 



D 



P^l + nksTss— (C + C* (a)) (a* - a*J + ... 



xy 



1 - 2 



C(a) 
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and 



Pyy = Pyy + ^^bTss [Ayy - 2a*-'P<^^^) (a* - a^) 



+ ... 



= P^'^+nkaT ( nn*~'^P*^'"'^ - 2a*~^ P*^'"''>\ ( n* - n* ) 

— I-yy -tUKBJ-SS y^a^g fyy j^^, ^^^-^ 1 ^ (^^^ ^^'^SS "sS^ 



= P^^ 

yy 



yy 



2g - C* ((^) _ 2 W — - 1 
C + iy* (a) + iC* (a) ) V a 



1-2 



v*{a) + C{a) \ (a 



C + y*{a) + \C*{(x)) 



^-1 +... 



and 



/ dT 

Pxx = DnksTss + (a* - a* J Dnks ( -q- 



{D-l)Pyy 

= DnksTss - 2 ( - 1 ) DnksTss - (-D - 1) Pyy + ... 
= DnkBTss-{D-l)P;i+2 

= Pf 



-I^n^.T. + (iP - 1) P- ( ^r^if Ifff, , 

V C + z/* (a) + |C (a) 



a 



- 1 



{C + y*{a) + \C*{a)){v*{a) + DC{a))\al, J 



The general moment equations are 



In the linear approximation, writing 



gives 



or 



which, for A/' > 2 is a simple linear equation for A}^ . 



(A13) 



(AM) 



(A15) 

(A16) 
(A17) 
(A18) 

(A19) 



(A20) 



(A21) 



(A22) 
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